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Abstract. We consider a recently reported micro-fluidic dye ring laser and study 
the full wave nature of TE modes in the cavity by means of finite-element simulations. 
The resonance wave-patterns of the cavity modes support a ray-tracing view and we 
are also able to explain the spectrum in terms of standing waves with a mode spacing 
5k = 2n/L c fi where L e g is the effective optical path length in the cavity. 
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1. Introduction 

Compact, efficient, and on-chip light-sources are of considerable interest for use in lab- 
on-a-chip applications [I] and recently there has been an increasing effort in realizing 
micro-ffuidic dye lasers based on glass or polymer j2j El IU El El Hj • 

Typically, the cavity designs rely on classical ray-tracing arguments rather than full 
wave simulations. In this paper we consider a geometry resembling that of Refs. El E] 
and offer a full wave study of the TE modes in the cavity. The resonance wave-patterns 
of the cavity modes support the ray-tracing view and we are also able to explain the 
mode-spacing of the spectrum in terms of standing waves. 

The paper is organized as follows: In Sec. 2 we present the geometry, in Sec. 3 
we address the mode spacing by quasi one- dimensional considerations, in Sec. 4 we 
numerically solve the wave equation, and in Sec. 5 we discuss aspects of optical gain. 
Finally, in Sec. 6 conclusions are given. 

2. Geometry 

We consider the two-dimensional laser resonator illustrated in Fig. ^ which corresponds 
to the planar cavities studied experimentally in Refs. |H1 El Ej. The cavity resembles 
a classical Fabry-Perot resonator and consists of two dielectric isosceles triangles with 
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Figure 1. Geometry consisting of a polymer-defined micro-cavity with an embedded 
micro-fluidic channel containing a dissolved laser dye. Out coupling of power from the 
cavity occurs through an evanescent-field coupling through an air gap to an adjacent 
polymer region where the output power is evaluated by an integral along the solid 
line 7. The dashed line indicates a typical optical path in the cavity. Simulations are 
carried out for point-source excitations at positions A, B, and C, respectively. 
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baseline £ po i and refractive index n pol separated by a microfluidic channel of width 4h 
containing a fluid with refractive index n c ^. Light is confined to the cavity by total- 
internal reflections at the polymer-air interfaces at an angle of incidence of 7r/4. Out 
coupling of power occurs through an evanescent-field coupling to an adjacent polymer 
region. In the experiments in Refs. [SHE] the microfluidic channel is filled by a dye doped 
liquid acting as gain medium. In Ref. i po \ ~ 700 /im, the cavity is pumped at the 
wavelength A = 532 nm by a pulsed frequency doubled Nd: YAG laser, and lasing occurs 
in the visible around A ~ 570 nm. For details on the pump power and lasing threshold 
we refer to Ref. 0. 

Throughout the rest of the paper we consider a typical structure with £ c h/£ P oi = 0.2 
and for the evanescent-field coupling we have w / £ po \ ~ 0.028 for the width w of the air 
gap. For the refractive indices we use n po i = 1.6 and n c h = 1.43. These numbers give 
an index step which is typical for a liquid and a polymer. However, we emphasize that 
the particular choice of numbers do not affect our overall findings and conclusions. 



3. Quasi one-dimensional approach to mode spacing 

We first estimate the mode spacing by considering a plane wave travelling around in the 
cavity, see Fig. ^ In this ray-tracing like approach we neglect reflections at the polymer- 
fluid interfaces which is justified by the very small Fresnel reflection probability 

R= ( n ^- n A ~Q,31% (1) 



,n po i + n ch; 

We imagine modes somewhat similar to whispering-gallery modes (WGMs) in resonators 
of circular shape. However, in this case the modes are subject to four total-internal 
reflections at an incidence angle of ir/A irrespectively of the mode-index and all modes 
have the same effective optical path length. Contrary to WGMs these modes have 
thus no cut-off for decreasing mode index caused by decreasing incidence angle. The 
accumulated phase during one round-trip of a plane-wave in the ring cavity is 

5(f) = kL cS + if (2) 
where k = 2tt/ X = uj/c is the free-space wave number, 

L c s = 2ri po i£ pol + 2n c h£ c h (3) 
is the effective optical path length in the cavity, and 



< <)S lfJ - \/Vr sin2 (f; 



<p = 4 x arg ( , = I (4) 

™s(f ) + y/ri£ - sin 2 (f ; 

is the phase picked up during the four total-internal reflections at the polymer-air 
interfaces at incidence angle of 7r/4. The resonance condition is 5<p = 2irm with the 
mode-index m being an integer. Obviously, the corresponding modes 
2ixm-ip 

k m = z (5) 

-^eff 
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Figure 2. Mode spectrum for a point-source excitation at position A, see Fig.^ The 
top panels show close-ups with a higher resolution of the respective peaks at m ~ 21 
and m ~ 27 indicated in the lower panel. 



are equally spaced with the mode-spacing 
2tt 
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4. Two-dimensional wave equation approach 

The full wave nature is governed by the wave equation [Hj 

V x V x E(r) = e(r)k 2 E(f) (7) 

where E is the electrical field and e(r) = n 2 (r) is the dielectric function. We solve 
the wave equation in a planar geometry for TE modes, i.e. E(f) = E z (f)e z and 
r = xe x + ye y . For the simulations we employ a finite-element method [Hj with "open" 
boundary conditions taken into account by perfectly matching layers (PMLs) at the 
edges of the simulation domain [TU] . see Fig. ^ This allows outgoing waves with 
negligible back reflection. 

We solve Eq. (JJJ) subject to a point-source excitation and modes are monitored 
by calculating the output power P on t(k) by integration along 7 in the polymer region 
adjacent to the cavity, see Fig. 1, for different values of k. The point-source has the 
appealing feature that it radiates isotropically in a homogeneous space and thus it will 
in general excite the full spectrum of cavity eigenmodes (except of course from the 
statistically few having a true node at the exact position of the point-source). 

In order to compare to the predicted mode spectrum, we have transformed the k 
values into a mode index 

m{k) = {kL cS + V )/2ti (8) 

and according to Eq. © we expect P out (m) to have resonances centred at integer values 
of m. Fig. 2 illustrates this in the case of a point-source excitation at point A, see Fig. 1. 

The over-all agreement between the full wave simulation and the quasi one- 
dimensional model is excellent, but from Fig. El it is also clear that the different peaks 
are slightly blue-shifted from integer values. The top panels illustrate this for two of the 
peaks indicated by green and red in the lower panel. This small shift may originate in a 
slightly modified phase shift at the edge with evanescent field coupling compared to the 
three other edges of the cavity. The small Fresnel reflection may also slightly modify 
the spectrum compared to the results derived from the quasi one-dimensional model. 

Figure El shows results in the range from m ~ 18 up to m ~ 29. When 
further increasing m the pattern of peaks persist with a small tendency that the 
peaks sharpen. This trend has been investigated up m ~ 100 where simulations turn 
highly computationally demanding (results not shown). However, since the quasi one- 
dimensional interpretation does not support a cut-off for increasing m we believe that 
a spectrum of equally spaced modes persist for increasing m. 

For decreasing m WGMs will typically experience a cut-off because the angles of 
incidence at some point do not support total-internal reflection. However, as discussed 
for the quasi one-dimensional model the particular class of modes in the present cavity do 
not share this property. In fact, in the simulations we have observed the modes down to 
m ~ 10 below which pronounced deviations from the quasi one-dimensional predictions 
start to emerge. Deviations most likely appear because the polymer-air interface has 
spatial variations on a length scale comparable to the wavelength of the light. In other 
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Figure 3. Electrical fields at m(fc) = 21.1000 and m(k) = 27.1355 for a point-source 
excitation at position A, see Fig. ^ 



words, the ray-tracing picture fails and concepts like total-internal reflection derived 
from Snell's law do not accurately capture the true wave physics. 

In order to verify that the peaks in Fig. |2] really do correspond to cavity modes we 
have studied the corresponding electrical fields at resonance, see Fig. El These fields 
resemble pure eigenfunctions of the resonator while off-resonance fields correspond to 
linear combinations of a larger number of eigenfunctions. Starting from e.g. the source 
point, the number of oscillations along one round trip equals m in full agreement with 
the quasi one-dimensional arguments. 

When the cavity is excited at different positions the overall output spectrum is 
the same such that peaks remain unshifted while changes are observed in the intensity 
distribution only. The reason is that different positions of the source will excite different 
linear combinations of eigenmodes (being correlated with the intensity level) while the 
eigenspectrum itself (being correlated with the resonance positions) remains unchanged. 
In Fig. 4 we illustrate this for different positions of the point source. The spectrum also 
reveals structure, though very broad with low intensity, in between integer values of 
m(k). This structure also corresponds to quasi eigenmodes which however are much 
more poorly confined to the cavity compared to the well-confined modes with integer 
values of m(k). 

5. Optical gain medium 

Lasing of course relies on the presence of an optical gain medium. In Refs. [3J|S] dissolved 
laser dye in the microfluidic channel provides the gain. While the dynamics of lasing 
is difficult to address we may with little effort investigate the influence of gain on the 
mode spectrum. The present numerical model allows for such investigations, but for low 
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Figure 4. Mode spectra for point-source excitation at positions A, B, and C, see 

pig. in 



concentrations general trends may be more easily analyzed with the aid of perturbation 
theory [TT] . Doping by e.g. Rhodamine 6G (Rh6G), as in the experiment [S], will change 
the refractive index in the channel accordingly, i.e. n<± + tiring where for the 

latter n Rh6G = n' Rh6G + m Rh6G . In the case of n Rh6G < n Rh6G < n ch we get 

k (E\Ae\E) 



Ak 



(^ch + ^Rh6c) 



n 



Rh6Gy 



(9) 



2 ( E \e\E) 

from which we expect a red-shift of the modes of the order Ak cx (^Rhec) 2 ^ a l° n g 
with a narrowing of the modes. A blue-shift may be observed in the case where 



n Rh6G > \/ n Rh6G( ?7 'Rh6G 



+ ^ch) 



6. Discussion and conclusion 



In this work we have used finite-element simulations to study the cavity mode spectrum 
of a micro-fluidic dye ring laser with a planar geometry resembling the one studied 
experimentally in Refs. |3] . We have performed a full wave study of the TE modes 
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in the cavity and found very good agreement with a quasi one-dimensional plane wave 
description with resonances corresponding to standing waves. 

In principle our simulations allow for an estimate of the quality factor of the 
modes, but realistic simulations for the experimental device require more details to 
be taken into account. For instance one would need to include the three-dimensional 
nature of the device to describe the radiation field accurately and the details of the 
evanescent field coupling would also influence the quality factor. Such issues add to the 
difficulty in addressing the dynamics of lasing so in this work we have only addressed 
the passive device. However, we have estimated a doping-induced shift of the spectrum 
by perturbative means. 

In the simulations we have considered mode-indices m(k) up to around 100 while 
in the experiments the corresponding typical mode index is estimated to be around 
two orders of magnitude larger. Nevertheless, we are confident that the standing-wave 
interpretation may be safely extrapolated to the experimental regime j3 13 E] due to 
the scale invariance of the wave equation [H] and the fact that this class of modes has 
no cut-off with respect to increasing mode index. 
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